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KNOT VISUALIZATION EXPERIMENTS FOR 
VERIFIABLE MOLECULAR MOVIES 



^q J. LI, T. J. PETERS, AND K. E. JORDAN 

u 

»-^H Abstract. Classical topological concepts are applied to under- 

^T stand high performance computing simulations of molecules writhing 

in three dimensional space. These simulations produce peta-bytes 
of floating point data, to describe 3 dimensional changes in molec- 
ular structure. A zero-th order analysis is achieved by viewing a 
computer animation synchronized with these changes. The per- 
formance demands for animation of this voluminous data can be- 
come problematic, but techniques from low-dimensional topology 
\[jJ are helpful. The 3D molecule is reduced to a lower dimensional 

model of a 1-manifold, which undergoes a piecewise linear approxi- 
mation for animation. An example is presented here to show how a 
C^ 1-manifold and its PL approximation could come to have different 

G embeddings as molecular writhing proceeds. This should serve as 

i *"^ | a cautionary warning to animators to respect established sufficient 

conditions for topological preservation so that the movies gener- 
,_^ ated will faithfully reflect the topology of the underlying model. To 

^. obtain this result, techniques from low-dimensional topology were 

,_^ joined with experimental mathematics and numerical analyses. 

in 

o 



1. Introduction 



For computations on protein molecules, a first conceptual reduction is 
provided by chemists through their notion of the backbone of the mole- 
cule (Figure [II . This is a skeletal structure which can be approximated 
by a 1-manifold, replacing a 3-dimensional model, to considerable compu- 
!— j tational advantage. This 1-manifold has a computational representation 

as a polynomial. For animation, this 1-manifold undergoes a piecewise 
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linear (PL) approximation. There are techniques to ensure that a static 
1-manifold and its PL approximation are ambient isotopic. 




Figure 1. The backbone of a molecule 



The primary example presented in Section [3] synthesizes results from 
two primary references [HI [T2] to advise caution for appropriate graph- 
ics approximations during dynamic visualization. It was motivated by 
two experiments conducted with our knot visualizing tool [12J. The first 
visual experiment is described in Section [2] to observe the convergence 
of a Bezier curve to its control polygon under collinear insertion, a new 
spline operation defined here. The second visual experiment is described 
in Section [3] where points were chosen and moved until the topological 
disparity appeared. 

The term 'molecular movies' has been coined to include ". . .cell and 
molecular animations ..." |18l ITS] . These animations can be created by 
individuals or can be programmatically driven. The shape preservation 
issues described here are applicable to both approaches. The creative 
animation process for individuals is supported by software |17| . which 
include shape changing techniques called 'morphing'. 

This preservation of topology is one property being explored under 
verifiable visualizations |9i: "That will consider both the errors of the in- 
dividual visualization component within the scientific pipeline and the in- 
teraction between and interpretation of the accumulated errors generated 
in the computational pipeline, including the visualization component." 

1.1. Related Work. There is contemporary interest HU E] to 
preserve topological characteristics such as homeomorphism and ambient 
isotopy between an initial geometric model and its approximation. The 
models here are 1-manifolds, represented by class of polynomials known 
as Bezier curves (Please see Section [2l , where these Bezier curves are de- 
fined by finitely many control points which become input to a standard 
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PL approximation algorithm known as subdivision. This subdivision al- 
gorithm is widely used in graphics and animation, so that a smooth curve 
is rendered from its PL control polygon. Sufficient conditions for a homc- 
omorphism between a Bezier curve and its control polygon have been 
studied |8j [13j [21] , while topological differences have also been shown 

The paper [12] presents techniques to create a class of examples - where 
a Bezier curve and its control polygon are not homeomorphic, as the 
Bezier curve is self-intersecting while the control polygon is simple, closed 
and equilateral, as shown in Figure [2] These were all static examples. 

The work [5] also established sufficient conditions for preservation of 
knot type during dynamic visualization of ongoing molecular simulations, 
as the knot type of the molecule is an important consideration for these 
biochemists. However, it relied upon a new approximation for each iter- 
ation of the perturbed curve. These approximation algorithms are more 
intensive than just continuing to perturb the initial PL approximation, as 
is examined here. 

Upper bounds were given |8j on the perturbations of the vertices or a 
broad class of PL approximations, so as to retain the same knot type of an 
underlying curve model representing the writhing molecule. It is easy to 
understand that drastic movements of the vertices of the PL model could 
create a different knot type in the graphics versus the perturbing curve, 
but it is shown that even small vertex changes beyond these established 
upper bounds could lead to incorrect knotting in the graphics display. Ap- 
preciating this subtlety is important to assess the pragmatic trade-offs of 
maintaining one PL approximation and perturbing it over multiple frames 
versus the more performance intensive approach of re-approximating the 
Bezier curve for each frame. 




(a) Simple C &c Self-intersecting C (b) Zoomed-in intersection 

Figure 2. Equilateral, simple C with self-intersecting C 
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1.2. Background. The strategy chosen is to perturb some control points, 
while creating a self-intersection of the Bezier curve, as described more 
fully in Section |1.3| An immediate difficulty encountered is illustrated in 



Figure [3J Figure 3(a) shows the Bezier curve from Figure [2] Figure [3 (b)| 
shows two resultant Bezier curves after a control point is perturbed to two 
different positions. Since both left and right images of Figure [3(b) | appear 
to be unknotted, it is possible that these unknots can be transformed into 
each other without ever passing through an intersection. This motivates 
us to create an example of perturbation changing the knot type of the 
underlying curve, which will be presented below. 




(a) Before Perturbations (b) Two Perturbations 

Figure 3. Perturbations 



1.3. Ideas. The key here is to prove the existence of self-intersection by 
demonstrating a change in knot type. Consider Figure [4] having the 
unknot indicated on the left, as the initial configuration. If the unknot 
is perturbed continuously to create the trefoil knot shown on the right 
of Figure |4l then a self-intersection must necessarily have occurred, as 
shown in the middle image of Figure HI The sufficient condition for the 
self-intersection is a change of isotopy class between the original and final 
knot. 

We generate a control polygon with a Bezier curve of high degree, by 
a new technique: collinear insertiorr] (described in detail in Section pi . 
Collinear insertion was developed by experimental mathematical methods, 
as assisted by a visualization tool, which is described elsewhere [12 . The 
experiments showed that perturbing one vertex of the subdivided control 

This terminology is used in other fields I16| . but with a very different meaning 
that should not cause confusion within this mathematical context. 
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FIGURE 4. Necessary self-intersection 



polygon typically left all other vertices fixed, so that embedding changes 
were not likely to occur. However, it was observed that moving the end 
point of a segment containing many collinear control points would cause 
more global changes. Collinear insertion was conceived to provide a close 
approximation between a Bezier curve and its control polygon (as happens 
after sufficient subdivision) while having the additional property of many 
collinear control points. 

With the carefully constructed example, we can perturb a vertex of the 
control polygon such that the control polygon remains simple during the 
perturbation, while the knotting structure of the Bezier curve changes (i.e. 
an unknot is perturbed to a nontrivial knot). Hence, the corresponding 
Bezier curve must become self-intersecting during this perturbation. The 
change of the knotting structure relies upon a numerical analysis, where 
the numerical tolerances can be easily adjusted to be within a single pixel 
resolution, merely by scaling the dimensions of the input according to the 
screen size, as clearly sufficient for computer animation. 

Admittedly, this is only a first order analysis but is likely sufficient 
for animation, but, more importantly, indicates the opportunity for more 
sophisticated numerical analyses, if deemed important for other appli- 
cations. For instance, the specific code implementations of the simplex 
method |10| and Horner's method [6 invoked in Section H could be sub- 
jected to a rigorous forward error analysis [23J to assess robustness, but 
such investigations are beyond the scope of this paper. Related error 
bounds on spline curves have appeared as interval splines [26 , which could 
afford alternative opportunities for other, more refined analyses. While 
the numerical details could vary, the topological perspective presented 
would adapt to those subtleties. 

1.4. Techniques: Collinear Insertion and Perturbations. Collinear 
insertion is similar to degree elevation [7] in that both methods produce 
higher degree Bezier curves. However, collinear insertion preserves the 
perimeter of the control polygon while changing the Bezier curve (see 
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Figures [5] & [6| , whereas degree elevation changes the control polygon but 
preserves the Bezier curve. 

Motivated by applications for dynamic visualization in a high perfor- 
mance computing environment, the preservation of topological integrity 
during perturbations has been investigated [3j |8]. For example, when 
visualizing a molecule twisting and writhing under local chemical and ki- 
netic changes, it is crucial that topological artifacts are not introduced 
by the visual approximations [8 . In particular, a perturbation on a con- 
trol point of a Bezier curve changes both the control polygon and the 
Bezier curve. The paper [3_ gives an upper bound on the perturbation of 
vertices of polyhedra to retain ambient isotopic equivalence. The paper 
|8„ Proposition 5.2] provides a sufficient condition to ensure that both a 
control polygon and its Bezier curve remain ambient isotopic during the 
perturbation. 

We present an example, consistent with previous upper bounds [3] on 
preservation of ambient isotopic equivalence, where a slight perturbation 
of a control point results in an ambient isotopic control polygon while 
changing the knot type of the Bezier curve. 

1.5. Outline. In Section[2] we describe the collinear insertion technique. 
In Section [3] we provide images of the simple, closed control polygons 
and Bezier curves constructed by the collinear insertion and perturbation. 
In Sections |4.1| and |4.2| we establish the knot types of the original and 
perturbed Bezier curves. In Section |4.3| we verify the non-self-intersection 
of the control polygon. 

2. Collinear Insertion for Bezier Curves 

Definition 2.1. [55] The parameterized Bezier curve, denoted as Bit), 
of degree n with control points P m e M. 3 is defined by 



B(*)= E ( n )t m (l-t) n - m P m ,t€[0,l], 



The PL curve determined by the points {Po,Pi, ■ ■ ■ , P n } is called the 
control polygon. If Pq = P„, then the polygon and the curve are closed. 
Otherwise, they are open. 

During the study of the existence of an equilateral control polygon 
with a self-intersecting Bezier curve, we tried to "slightly" perturb a ver- 
tex so that the control polygon would remain simple, but simultaneously 
the Bezier curves would differ "dramatically" by a change of knot type. 
However, for many examples, if a vertex is "slightly" perturbed, then the 
Bezier curve is "slightly" changed, i.e. the knot type remains the same. 
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One may use the de Casteljau algorithm to create a new control poly- 
gon which is closer [20 to the Bezier curve in Hausdorff distance, but 
this also shortens the length of edge of the control polygon, such that the 
perturbation of a vertex only effects a small region of the control poly- 
gon. Consequently, the perturbation may not change the Bezier curve 
sufficiently to change the knot type. 




FIGURE 5. Successive Collinear Insertions 

A fundamental intuition about Bezier curves is that the curve tends 
towards the control points. It follows that adding more collinear control 
points should draw the curve closer to those edges. We were able to cor- 
roborate that visually through experiments, as shown in Figure [5] where 
the sequence of smooth curves are obtained from successive collinear in- 
sertions. The innermost curve is defined by four control points, whereas 
the curve closest to the control polygon has 25 control points, with the 
intermediate ones defined by 7 and 13 control points, respectively. That 
experimental evidence was sufficient for creation of this example, while 
leaving a formal convergence proof as a subject for further investigation. 

The illustrative example of Figure [5] is of an open curve. The rest of 
this section treats closed curves. Suppose that there is an initial set of 
control points, denoted as {vq, v±, . . . , i> n -i, ^o}i defining a closed curve of 
degree n. (Note that the vertices {w , i>i, . . . , w n -i} with v ^= v n -i would 
define an 'open' curve of degree n — 1. ) Now we describe the procedure 
to raise the degree by adding collinear control points. We designate this 
as collinear insertion. 

(1) Take the midpointaj along each line segment. 



For ease of exposition, the collinear points inserted are all midpoints. However, it 
is clear that other schemes for insertion of collinear points are possible. This is similar 
to the typical exposition of the de Castlejau algorithm |7|, where the midpoints are 
used to explain the concept, but other intermediate points are equally possible. 
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(2) Add these midpoints to the initial set of control points to form a 
new set of control points: 

r w + «i v n -i + v 
{vo, — 2 — ,«!,..., W„_i, ,v }. 

Consequently, the perimeter of the control polygon remains the 
same, but the Bezier curve changes. (Retaining this original 
perimeter, even as more control points are added, supports an 
easy method to perturb many control points simultaneously while 
only perturbing the original vertices.) The degree is raised from 
n to 2 * n. 

(3) If the process is repeated k times, then we obtain a set of 2 fe * n 
control points, which determines a closed Bezier curve of degree 
2 k *n. 




(a) Initial curves (b) Higher degree curves 

Figure 6. Raising the degree 



Figure 6(a) shows a closed control polygon with a closed Bezier curve 
of degree 7 determined by 7 vertices, and Figure [6(b)| shows the control 
polygon with a Bezier curve of degree 112 determined by 112 control 
points, constructed by the collinear insertion, which is the focus of the 
rest of this paper. 

3. The Defining Data for the Example 

We start with 7 control pointf|j listed below to determine a closed 
control polygon. 

(1.9817, -1.7646, -4.5897), (-1.3841185, 4.6825505, 0.913541), 

At least four decimal digits are provided here. If higher accuracy is desired, more 
digits may be chosen. 
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(-3.2983075, -4.0566825, 2.686189), (-0.1232995, 2.768254, -2.463584), 

(3.9079915, -4.533357, 1.2263705), (-3.935983, -0.438272, -0.983365), 

(3.218174, 4.296123, 2.1124595). 

We then use the collinear insertion to insert midpoints. Repeating 
this process four times, we obtain 112 control points. The corresponding 
control polygon and Bezier curve of degree 112 are shown in Figure [7(a) 



which repeats Figure |6(b)| for a visual comparison to Figure |7(b)| which 
is obtained by perturbing one of the original vertices, as described below. 
Denote the control polygon and the Bezier curve in Figure 7(a)| by V\ and 




(a) Initial polygon and curve (b) Perturbed polygon and curve 

Figure 7. Perturb v to v' 

Denote the vertices of V\ (the original control points) as vq, Vi, . . . ,Vq. 
Note that the vertices adjacent to vq are v\ and vq. We fix v\ and vq 
while moving the line segments VqV\ and VqVq, by perturbing vq from 
(1.9817, -1.7646, -4.589 7) to (1.3076,-3.3320,-2.5072), denoted by v' , 
as shown in Figure [7(b) | Denote the resultant control polygon and Bezier 
curve as Vi and Z?2- 

The perturbation of vo to v' a is an affine transform, performed by the 
elementary graphics operation of rotation of a point about an arbitrary 
line [25]. The line of rotation is formed by the points v\ and vq. This 
creates a new control polygon, where only two segments have been altered, 
those containing v' , which are Viv' Q and Vqv' . 

When vq is perturbed to v' , the control polygon remains simple, but 



the unknot B\ is changed to a nontrivial knot B2 (See Section 4.2 1. The 



different knotting structures between B\ and B2 imply that there must 
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be an instant during the perturbation, when the Bezier curve is self- 
intersecting. This is an example of a perturbation which preserves the 
ambient isotopy of the control polygon but changes the knot type of the 
Bezier curve. 

4. Topological Differences 

The following three Subsections show that 

• the initial Bezier curve is unknotted. 

• the perturbed Bezier curve is a nontrivial knot, indicating that 
a self-intersection of the Bezier curve must have occurred during 
this perturbation, and 

• the control polygon remains simple throughout the perturbation. 

We prove the knottedness and unknottedness under the sub-pixel cri- 
terion, that is, we accept results generated by numerical methods if any 
numerical tolerances can be adjusted to provide sub-pixel resolution, as 
sufficient for this context. 




(a) Initial curve (b) Perturbed curve 

Figure 8. UnKnot and Knot 



4.1. Initial Unknottedness. For the unknottedness of Si, consider the 
projection shown in Figure 8(a)| obtained by taking each z-coordinate to 
be zero, that is, the projection onto the xy-plane. Denote the projected 
Bezier curve by C2d(t). Define a function fnS(t\,t 2 ) ■ [0, 1] x [0. 1] — >• R 
by 

fnS(h,t 2 ) = \\C2d(ti)-C2d(t 2 )\\. 

If fnS(ti,te) = and t\ ^ t 2 , then we have C2d(t 1 ) = C2d{t 2 ) and t\ ^ 
t 2 , which implies that the projected Bezier curve is self-intersecting, and 
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the self-intersection has t\ and t 2 as parameters. In other words, we seek 
the roots of fnS(ti,t2) with t\ ^ £2 to determine the self-intersections. 
Note that fnS{ti,t<}) is a polynomial of degree 112. We use the simplex 
search method [10 implemented by Mathlab function 'fminsearch' to find 
the roots of this high degree polynomial. 

The pairs of parameters where the self-intersections occur are listed 
below, which correspond to intersections labeled as '1, 2, 



Figure 8(a) 



[0.0488, 0.4614], [0.0861, 0.7918], [0.3473, 0.6931], [0.5126, 0.9915]. 
Using these par ameter values to evaluatajthe original Bezier curve in 



parameter values: 



3D (Definition 2.1 1, we get the pairs of points corresponding to the above 



[(0.8309, 0.4397, -2.7081), (0.8308, 0.4397, -1.2072)]; 

[(-0.0435, 2.0929, -1.2807), (-0.0435, 2.0928, 0.6747)]; 

[(-1.8672, -1.0031, 0.4288), (-1.8672, -1.0032, -0.3359)]; 

[(2.0548, -1.4062, -0.3480), (2.0548, -1.4061, -4.1932)]. 



Note that the orientation of B\ is shown by arrows in Figure 8(a) By 
comparing the z-coordinates of points listed above, we find that '1' and 
'2' are under crossings, while '3' and '4' are over crossings. It follows that 
B\ is simple, within the numerical limits imposed by this finite precision 
data and these floating point computations. 

Observe specifically that the sub-curve from '3' to '4' is over the other 
portions. So using the published notation [15], a Reidemeister move of 
Type 2b will eliminate the crossings '1,3' and '4', by pulling the sub-curve 
from '3' to '4' downwards. Consequently the resultant diagram after the 
Reidemeister move has only one crossing, which implies B\ is unknotted. 

4.2. Nontrivial Knot Occurrence. For the nontrivial knottedness of 
B2, we start from a projection obtained by taking each z-coordinate to be 
zero, as shown in Figure [8 (b)| 

Again, using the MathLab function 'fminsearch', we find the following 
pairs of parameters where the projected curve is self-intersecting: 

[0.0761,0.4240], [0.0928,0.7830], [0.3473,0.6931], [0.5039,0.9575]. 
Evaluating (aga in, using Horner's method) the original Bezier curve in 



3D (Definition 2.1) at these parameters produces the corresponding pairs 
of points: 

[(-0.1265, 0.9308, -0.6850), (-0.1266, 0.9308, -1.2854)]; 



This evaluation is performed by Matlab, where Horner's method [5] is invoked. 
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[(-0.4389, 1.8160, -0.2901), (-0.4389, 1.8159, 0.5159)]; 

[(-1.8672, -1.0032, 0.4289), (-1.8671, -1.0032, -0.3358)]; 

[(1.8761, -1.0622, -0.5163), (1.8761, -1.0623, -1.1327)]. 

The under and over crossings follow from comparing the z-coordinates 
in each pair, as shown by the diagram in Figure [9] 




Figure 9. Diagram 



In the diagram, the crossings are labeled as Ci, . . . , C4 and regions are 
ro, . . . ,r$. According to [T], we obtain the following corresponding matrix 
of the diagram: 

" 
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-1 


-t 
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-t 
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t 





-1 





1 







1 





-1 


t 



t 

Deleting any two consecutive columns of the above matrix, computing 
the determinant and normalizing the polynomial, we get the Alexander 
polynomial 1 — 3i + t 2 , which implies the knot is nontrivial [1] . 

4.3. Control Polygon Remains Simple. Let a be the angle of rotation 
from vq to v' . Let v{9) be the point between vq and v' during the rotation, 
where 9 e [0, a] and it indicates the angle of rotation from vq to v(9). This 
rotation function is the standard rotation matrix from an introductory 
graphics course [25 . Now, to show the non-self-intersection of the control 
polygon during the perturbation, it suffices to show that the line segments 
v(9)v\ and v(Q)vq do not intersect with other line segments of the control 
polygon, for all 9 £ [0, a]. 
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For each i = 1, • • • ,5, consider v(6)vi and ViVi+i- Parametrize the line 
segments by: 



v{0)v x : v(8) + ( Vl - v(0))s, s e [0, 1] 
ViV i+1 : Vi + (v i+ i - Vi)t, t <E [0, 1] 

And then consider the equation 

v(9) + (vi - v(6))s = v % + (vi+i - Vi)t. 

It has no solution for 9, s and t with 9 € [0, a] and s, t € [0, 1], so v(9)vi 
does not intersect ViVi+i. 

We can similarly verify that v(9)v$ does not intersect ViVi+i for each 
i = l,--- ,5. 

5. Conclusion 

In this work, classical topology is complemented by experimental math- 
ematics and numerical analysis to provide insight into pathologies that 
could occur when creating synchronous visualizations of molecular simu- 
lations. A representative example is presented where a small data pertur- 
bation produces different knot types between a polynomial curve and its 
PL approximation used for graphics rendering. The disparity exhibited 
serves as a cautionary note to the application community for rigorous 
attention to graphics approximations during dynamic visualization. The 
analyses assume a sub-pixel criterion, but, more importantly, the methods 
(including collinear insertion) are likely to serve as a catalyst for domain 
scientists to investigate the relationship between smooth and PL knotting 
structures. 
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